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ABSTRACT 

Context. Recognizing the transport of Energetic Neutral Atoms (ENA), from their place of birth to Earth orbit, has become an 
important issue in light of the forthcoming launch of the NASA SMEX mission IBEX, which is devoted to imaging of the heliospheric 
interface by in-situ detection of ENAs. 

Aims. We investigate the modifications of both energy of survival probability of the hydrogen ENA (H ENA) detectable by IBEX 
(0.01 - 6 keV), between the termination shock and Earth orbit. We take into account the influence of the variable and anisotropic solar 
wind and of solar EUV radiation. 

Methods. Energy changes of the atoms are calculated by numerical simulations of the orbits of H ENA between ~ 100 AU from the 
Sun and Earth orbit, taking into account solar gravity and Lyman-a radiation pressure, which is variable in time and depends on the 
radial velocity of the atom. To calculate the survival probabilities of the atoms against ionization, a detailed observation-based 3D and 
time-dependent model of H ENA ionization is constructed, and with the use of this model the probabilities of survival of the atoms 
are calculated by numerical integration along the previously-calculated orbits. 

Results. Due to radiation pressure, H ENA reach the Earth orbit practically without energy and direction change, apart from the atoms 
of energy lower than 0.1 keV, during high solar activity. The survival probability of H ENA increases from just ~ 2% for the slowest 
detectable ENA at solar minimum to ~ 80% for the fastest ENA. For a given energy at Earth orbit we expect fluctuations in the survival 
probability of amplitude between ~ 20 percent at 0.01 keV to just a few percent at 6 keV and a modulation of survival probability as 
a function of the location at Earth orbit, ecliptic latitude of the arrival direction, and phase of solar cycle with an amplitude of a few 
dozen percent for 0.1 keV atoms at solar minimum to a few percent for 6 keV atoms at solar maximum. 

Conclusions. With appropriate account of local transport effects IBEX should be able to discover departures from symmetry in the 
flux of H ENA from the heliospheric interface at a level of a few percent. 

Key words. Interplanetary medium - ISM: atoms - (Sun:) solar wind - Sun: UV radiation 
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1. Introduction 

Hydrogen Energetic Neutral Atoms (H ENA) provide impor- 
tant information about the physical state of plasma in re- 
mote part s of the heliosphere, especially in the heliospheric 



, ospnenc 

interfac e (lGruntmanlll992l: iH sieh et all 1 1992!: lGruntmanlll997t 



Gruntm anet all r20()U" [Sch erer & Fahr 2003: iFahr & Scherei 
2004: Heerikhuis en et alj l2007i: fahr et al.i i2007l: ISternal et al. 
2008). Measurements of heliosph eric H ENA are planned by 
the NASA SMEX mission IBEX dMcComas et al.ll2004 12005L 
12006.) . IBEX will be a spin-stabilized Earth satellite observing 
Energetic Neutral Atoms (ENA) in the 0.01 - 6 keV energy 
band, with sensors looking perpendicularly to the spin axis di- 
rected at the Sun. The prime target of the mission will be imag- 
ing of the heliospheric interface via neutral atoms and the goal 
of an individual observation will be to obtain the kinematic pa- 
rameters of the locally-registered atom beyond the termination 
shock (TS). To that end, understanding of the processes affect- 
ing the atoms between their source region and the detector is 
needed. In particular, a control on the modifications of the kine- 
matic parameters of the atoms between the source region and 
the detector and a grasp on loss processes along the way are 
necessary. So me insights into this problem were presented by 
iBzowski & Tarnopolskl(2006.) . 
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Since the ENA gas inside the TS is collisionless at all en- 
ergies, then - from Boltzmann equation - the transport of ENA 
from TS to the inner heliosphere is governed solely by dynamic 
effects and gain and loss processes. The dynamic effects include 
heliocentric acceleration by a combined action of the solar grav- 
ity and the opposing solar radiation pressure, whose strength de- 
pends on the radial velocity of the incoming atom. The losses 
include all known ionization processes operating in the helio- 
sphere: charge exchan ge, EUV ionization, and electron impact 
(iRuciiiski et alj [19961) . For inwardly-traveling ENA, the gains 
(considered by IBEX as unwanted foreground) include charge 
exchange between CME, CIR, and pickup ions on the one hand, 
and all kinds of neutral populations present inside TS on the 
other hand, as well as production of ENA in the plasma environ- 
ment of the spacecraft. In this paper, they are neglecte d because 
they will constitute a small addition to the net signal dFahr et al.l 
2007). 

The paper discusses modifications of kinematic parameters 
and survival probabilities of H ENA, suitable for detection by 
IBEX, traveling from the heliospheric termination shock (TS) 
to 1 AU. First, we discuss the radiation-pressure force affect- 
ing the atoms and the resulting trajectories of individual H ENA 
reaching perihelia of their trajectories at Earth orbit. Then we 
briefly present the ionization processes destroying the incom- 
ing H ENA will be briefly presented, providing a longer dis- 
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cussion in the Appendix. Subsequently, we show the survival 
probabilities of the H ENA, with emphasis on asymmetries and 
anisotropies that result from local conditions in the solar wind. 
The discussion is based, as much as possible, on available solar- 
wind and solar-radiation output data and does not consider mea- 
surement technicalities, such as detection efficiencies, geometric 
factors, proper motion of the spacecraft etc., which are left to a 
future study. Emphasis is placed on atoms at the energies ex- 
pected in the heliosheath region, i.e. up to 0.2 keV; faster atoms 
are not, however, disregarded, since the studies mentioned ear- 
lier predict a non-vanishing production of H ENA to the upper 
limit of the IBEX sensitivity threshold. Attention is drawn to 
solar-cycle-related effects, which include changes in the solar- 
wind anisotropy and net flux of solar Lyman-a pressure, and to 
fluctuations on a timescale of the solar rotation period. 

The term "source region" frequently used in the paper should 
be understood as an origin of the H ENA, sufficiently distant 
from the Sun for its dynamical influence on the atoms to be neg- 
ligible. It was arbitrarily set to be ~ 100 AU from the Sun. An 
atom that is unaffected by any radiation pressure and at infinity 
travels at 12 km/s, at 1 AU will have accelerated to 43.8 km/s, 
while an atom that travels with the velocity 12 km/s at 100 AU 
from the Sun, at 1 AU will have accelerated to 43.6 km/s and 
the difference between the velocities is then just 0.2 km/s; when 
radiation pressure is added and larger velocities are adopted, this 
difference will decrease further. 

Dynamical conclusions discussed in Sect. 2 are valid for all 
atoms reaching perihelion at 1 AU from the Sun irrespective of 
exact location in the 3D space; the conclusions about the sur- 
vival probability, however, presented in Sect. 3 were obtained 
for H ENA reaching perihelia of their trajectories at 1 AU in the 
ecliptic plane, which is referred to as "at Earth orbit". 

2. Dynamics of H ENA atoms between the 
termination shock and inner heliosphere 

2.1. Radiation pressure and equation of motion 

The velocities of heliospheric H ENA that can be observed by 
IBEX are in the approximate range ~ 50 - 1000 km/s, which 
corre s ponds to an energy range 0.01 - 6 keV (McComas et al.l 
12004 l2005Ll2006h . The motion of such atoms between the ter- 
mination shock and inner heliosphere is determined on the one 
hand by the solar gravity force, and on the other hand by res- 
onance interaction with the solar Lyman-ff radiation. The latter 
phenomenon produces a force called radiation pressure, which 
decreases with heliocentric distance following 1/r^, similar to 
solar gravity and has a radial direction that is away from the 
Sun. Hence, it is convenient to express the radiation pressure 
as a factor fi of (over)compensation of solar gravity. The solar 
Lyman-a line is about 0. 15 n m wide and features a self-reversed 
profil e with a central trough jLemaire et alj|2002l ; iLemaire et alj 
12005 ). Given the non-flat shape of the solar Lyman-a line (see 
Fig. [TJ on one hand, and the range of radial velocities of the in- 
coming H ENA atoms on the other hand one expects that, due 
to the Doppler effect, the radiation pressure acting on the atoms 
will be a non-negligible function of their radial velocities with 
respect to the Sun. If the heliospheric gas is optically thin, all 
neutral atoms in the heliosphere obey the equation of motion in 
the form: 



Table 1. Parameters of the model of radiation-pressure depen- 
dence on radial velocity Vr, expressed in km/s, and on total flux 
/tot, expressed in cm"^ s defined in Eq.©. 



d^r/dfi ^-GM[l-^ (v, (r) , /tot (f, 0))] r/|rp 



(1) 



A ■■ 

D 

H 



■■ 2.4543 • lO-'', 
= 0.73879, 
= 0.47817, 



B ■■ 

F 

P 



: 4.5694 
■- 4.0396 
: 4.6841 



10-^ 
10-2, 



C 
G 
Q 



3.8312 • \Q-\ 
3.5135 • 
■■ 3.3373 • 10-" 



the product of the gravity constant and the solar mass, /tot (f, 4>) is 
the wavelength-integrated solar Lyman-a flux at a time t and he- 
liolatitude (p, and // is the radiation pressure exp ressed as a com- 
pensation factor of sol ar gra vity. According to IScherer & Fahri 
(1996) and Ouemerais (2000), the optical depth of the neutral in- 
terstellar gas inside the heliosphere cannot be neglected beyond 
~ 10 AU from the Sun in the case of atoms of thermal ener- 
gies (<~ 0.01 keV); for radial velocities exceeding ~ 30 km/s, 
however, the gas is optically thin and the above-mentioned ap- 
proximation is valid. 

Based on a fit of Eg. Q to 9 sola r Lyman-o- line pro- 
files observed by iLemaire et alj ( l2002h . the solar Lyman-ff 
line profile can be parameterized by a simple function of 
wa velength- and disk-integrated flu x /tot (Tarnopolski, 200 8, the- 
sis; [^I^ioisE^liowsy (|2003); lio^^ 

H (vv, /tot) = A (1 + B /tot) exp (-Cv^) (2) 
X [ 1 H- D exp {-Fvr - Gvfj + H exp [Pv,. - gv^)] 

where A,B,C,D,F,G,H,P, and Q are parameters of the fit, 
which are compiled in Table [T] Conclusions presented in the 
following part of the paper are based on numerical solutions of 
the equation of motion priovided in Eq. ([T) in connection with 
Eq. (O, which were obtained with initial conditions including 
velocity vectors perpendicular to the local radial direction at a 
starting point located at Earth orbit, with velocity magnitudes 
corresponding to energy channels planned for the IBEX detec- 
tors. 

Observations of /tot have been performed for about three 
sola r cycles and inev itable gaps were filled with proxy val- 
ues (iTobiska et al.l !2000). The absolute calibration is still some- 
what uncertain. In pre sent paper, as basi s for /tot the SOLAR 
2000 model was used JWoods et al.ll2000h . which was also em- 
ployed by'Lemair e et al.l (l2002h to calibrate their observed solar 
Lyman-Qf line profiles. In the calculations the Lyman-o' flux was 
assumed to be constant in time and for solar minimum to be 



where r (f) is the position vector of the atom with respect to the 
Sun at a time f, - (dr/dt)-r/\r\ is its radial velocity, and G M is 



equal to /min - 3.7 x lO" cm ^ s ' and for solar maximum to 
7,^^^ = 6.0-10" cm-2 s-i. 

2.2. Trajectories of H ENA and modifications of their 
energies between fhe source region and 1 AU 

Because of the geometry of its lines of sight, IBEX will measure 
only the inward-traveling atoms that reach perihelion at 1 AU. 
Such atoms are detected at approximately zero radial velocity 
with respect to the Sun, but for most of their paths they have ra- 
dial velocities that are almost equal to their net velocities. If their 
net velocity in the source region is larger than the spectral range 
of the solar Lyman-a line, the atoms do not experience any ra- 
diation pressure underway, aprat from within a few AU from the 
Sun, where, before reaching perihelion, they their radial veloci- 
ties will be reduced to zero (see Fig.|2]). Since the velocities are 
large, however, the trajectories will be weakly modified by so- 
lar gravity and remain close to straight lines. On the other hand, 
slower ENA (whose velocities remain within the spectral range 
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Radiation pressure for H ENA Solar hy-a Une-integrated flux: 



monthly averages 




-200 -100 100 200 1990 1992.5 1995 1997.5 2000 2002.5 2005 

radial velocity (km/s) time (yr) 

Fig. 1. Radiation pressure acting on H atoms, expressed as a factor fi of compensation of the solar gravity force, shown as a function 
of their radial velocity with respect to the Sun for solar minimum and maximum (left-hand panel), and wavelength- and disk- 
integrated monthly-averaged solar Lyman-a flux based on the SOLAR 2000 model (right-hand panel). The profiles of radiation- 
pressure factor fi are shown according to the model defined in Eq.(|2]i for the net flux = 3.7 ■ 10" cm"^ s ' (solar min.) and 
■^max = 6.0 ■ 10" cm~^ s~^ (solar max.). 



Radiation pressure before perihelion (solar min.) 



Radiation pressure before perihelion (solar max.) 
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Fig. 2. Instantaneous radiation pressure acting on H ENA during 81 days before perihelion located at 1 AU, expressed as the 
gravity-compensation factor fi and shown for a few energies corresponding to the lowermost IBEX energy channels and equal to 
0.01 keV (solid), 0.015 keV (dots), 0.029 keV (dash-dot), 0.056 keV (dashed), 0.107 keV (dash-dot-dot), 0.208 keV (lowermost 
solid lines), for solar minimum (left-hand panel), and maximum conditions (right-hand panel). The perihelion velocities correspond, 
respectively, to 43.8, 53.6, 74.5, 103.6, 143.2, and 199.6 km/s. 



of the solar line) will be stronger affected by radiation pressure. 
Effectively, almost all changes in the kinetic energy of incoming 
atoms occur within ~ 3 solar rotations before the 1 AU peri- 
helion or an even shorter time (the faster the atom moves, the 
shorter is the effective interaction time), as illustrated in Fig. [3] 
During solar minimum, the effective radiation pressure is, how- 
ever, close to /i = 1, which produces little modification of the 
trajectories and energies. In contrast, the trajectories of the slow- 
est atoms in the discussed energy range will be affected con- 
siderably during solar maximum, when radiation pressure over- 
compensates significantly for solar gravity. Modification of their 
energies between the source region and perihelion at 1 AU is 



shown in Fig. |4] and in the left-hand panel of Fig. |5] The fig- 
ures demonstrate that energy modifications are negligible for all 
IBEX energy intervals apart from the slowest atoms during solar- 
maximum conditions: in this case, they can reach 50%. Atoms 
faster than 200 km/s at 1 AU practically do not change their 
energy. The difference between the arrival direction and the di- 
rection at the source region is negligible at solar minimum (see 
Fig.©, but at the lowest energies increases to 14° at solar maxi- 
mum. Deflections larger than the 7° of the IBEX pixel width are 
therefore possible only in the case of energies below ~ 0.03 keV. 



4 



M.Bzowski: Modifications of H ENA between termination shock and Earth orbit 



Change of kinetic energy before periheUon, Change of kinetic energy before perihelion, 

solar min solar max 
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Fig. 3. Ratios of energy of H ENA at time f before perihelion located at 1 AU to the energy in the source region during 8 1 days 
before perihelion for a few 1 AU energies, corresponding to the lowermost IBEX energy channels and equal to 0.01 keV (solid), 
0.015 keV (dots), 0.029 keV (dash-dot), 0.056 keV (dashed), 0.107 keV (dash-dot-dot), 0.208 keV (sohd hnes), for solar minimum 
(left-hand panel) and maximum conditions (right-hand panel). 



Source/1 AU energy ratio for solar min. and max., 
equatorial and polar directions 



Pole/equator source energy ratios, solar min. and max. 
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Fig. 5. Left-hand panel: Ratios of H ENA energy in the source region fis.eq to their energy at 1 AU for atoms incoming from 
polar regions Es/Eia\j (broken lines) and from equatorial regions /ieq/^iAU (solid lines) during solar minimum and maximum 
conditions. The upper pair of lines corresponds to solar maximum, the lower pair to solar-minimum conditions. Right-hand panel: 
ratios of energies in the source region of atoms incoming from the poles and from an equatorial band Es/Eeq, shown as a function 
of their energy at 1 AU for solar minimum and maximum conditions (respectively, solid and broken line). 



As is evident from Fig. \T\ the solar Lyman-o- flux ex- 
hibits appreciable fluctuations on both monthly and daily scales. 
Analysis of daily time series, split into Carrington-rotation in- 
tervals, demonstrates that the ratios of standard deviations of the 
Carrington-averaged values to the Carrington-averaged values 
vary from ~ 2% at solar minimum to ~ 5% at solar maximum. 
Based on the results presented above, it is surmised that such 
fluctuations should not introduce appreciable fluctuations in the 
energies and directions of the incoming ENA atoms. 



2.3. Effects of latitudinal anisotropy of the Lyman-a flux 

Potentially stronger variation than the fluctuations discussed in 
the previous section might be expected because of the lati- 
tudinal anisotr opy in the solar Lyman-ff output. As proposed 
by Cook et all (Il980i Il98ll) and observed directly by Auched 
(2005) during solar minimum and indirectly by Prvor et aT| 
(119921) during solar maximum, the line- and disk-integrated solar 
Lyman-of flux is a function of heliolatitude. The theory and lim- 
ited observations available appear to imply that the polar flux /poi 
should be equal to about 0.8 of the equatorial flux /equ and that 
this ratio should not change appreciably during the solar cycle. 
More research is certainly required to establish details of this 
phenomenon and the anisotropy-related conclusions presented 



M.Bzowski: Modifications of H ENA between termination shock and Earth orbit 



5 



Energy at 1 AU vs energy at TS 
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Fig. 4. Energy of H ENA at 1 AU versus energy in the source 
region at 100 AU. Solid line: solar-minimum conditions, bro- 
ken line: solar maximum-conditions, dotted line: purely bal- 
listic case with no radiation pressure, shown for comparison. 
Projections of dots on the horizontal axis mark the IBEX en- 
ergy channels, apart from the lowermost one, which marks the 
low-energy sensitivity limit. 



was modified to be 



yu (0) = yu (0) y cos2 ^ + p sin^ (3) 

Simulations demonstrated that for a selected H ENA energy at 
1 AU (denoted /siau), the modulation of the ratio of H ENA 
energies, in the source region at the ecliptic poles Ef^/E^, as a 
function of ecliptic longitude of perihelion for the lowest ener- 
gies detectable by IBEX Eipjj - 0.01 keV, is approximately 
4%, and at 0.06 keV, is already as small as 1%. This effect can 
therefore be considered to be negligible. 

Irrespective of the solar-cycle phase, the source-region en- 
ergies iisi-c (0) of the H ENA arriving at 1 AU from differ- 
ent latitudes (f> with the same energy Ei au, vary with (p by no 
more than ~ 10% in the case of lowest energies (see the right- 
hand panel of Fig. |5]l. These variations decrease to below 5% at 
E\ AU - 0.02 keV and disappear altogether for £1 au > 0.2 keV. 

The eff'ects of anisotropy in the solar Lyman-a flux persist 
both during solar minimum and solar maximum. During solar 
minimum, the atoms originating in the polar regions are trans- 
mitted with a very small energy change for the entire energy 
range (see the broken line in the solar-minimum pair of lines 
in the left-hand panel of Fig. |5]l, while a 20% change in the 
energy of the slowest atoms originating in ecliptic latitudes is 
expected (the solid line). Similarly during solar maximum, the 
energy change for polar directions is smaller than for equatorial 
directions (see the upper pair of lines in the left-hand panel of 
Fig.E). 



Deflection at 1 AU from direction at source 
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Fig. 6. Difference between the original direction of ENA within 
the source region and the direction at 1 AU for solar minimum 
conditions (solid line) and solar maximum (broken line) as a 
function of energy at 1 AU. 



below must be treated with caution. By adopting this anisotropy 
and assuming that the symmetry plane of the latitudinal pattern 
of the Lyman-a flux is the solar equator and the line profile does 
not change with latitude, we would expect a yearly modulation 
of the flux of H ENA originating in polar regions (because of the 
Earth travel in heliolatitude by +7.25°) and a differentiation of 
H ENA energies as a function of arrival-direction helioaltitude. 

To simulate the "flattening" of the net Lyman-or flux at polar 
latitudes, the radiation-pressure compensation-factor fi in Eq.(|2|i 



2.4. Travel times of H ENA from source region to 1 AU 

For the high end of the IBEX range of H ENA energies at 1 AU, 
the travel times from the source region are practically equal to 
the travel time in the force-free case (straight-line, constant- 
speed motion). The differences are approximately 1 solar rota- 
tion period in duration for energies between 0. 1 and 0.2 ke V, and 
for more energetic atoms are even smaller For ~ 5 keV atoms, 
the travel time is equal to 6 Carrington rotations and for 0.2 keV 
atoms, it increases to 32.5 rotations. Deviations from the force- 
free case begin for the 1 AU energies about 0.05 ke V and become 
dramatic towards lower energies, as shown in Fig.|7] An excep- 
tion is the purely-ballistic case (which is defined here with no 
radiation pressure); noticeable differences can be seen already at 
0.1 keV of El AU energy and become enormous towards lower 
energies in the IBEX sensitivity band, for which the travel times 
of such atoms are twice as long as those of H ENA. The differ- 
ences in travel time of H ENA due to different values of net solar 
Lyman-ff flux during the solar cycle, range from approximately 
10% at 0.06 keV to 20% at 0.01 keV. At solar maximum, when 
the net flux is higher, the atoms detected at a given energy require 
less time to arrive at 1 AU than during solar minimum, although 
the radiation pressure is higher. This paradox is explained by the 
fact that both radiation pressure and solar gravity are effective 
only relatively close to the Sun. During solar maximum condi- 
tions, the atoms become decelerated just prior to their arrival at 
1 AU, and consequently spend most of their travel moving much 
faster than the atoms that arrive during solar minimum, when 
the atoms follow almost force-free trajectories. This effect, how- 
ever, exists only for the lowest-energy atoms below ~ 0.06 keV 
at 1 AU. 

If the solar Lyman-or flux is indeed lower at the solar poles 
than at the equator, as discussed in the preceding section, we 
would expect differences in the travel time between the equa- 
torial and polar regions, which for the lower end of the IBEX 
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Travel time from 100 to 1 AU, solar min. and max. 
equatorial and polar trajectories 
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Fig. 7. Left-hand panel: Travel times of H ENA from the source region at 100 AU from the Sun to 1 AU as a function of energy 
at 1 AU. The uppermost solid line corresponds to the travel time with no radiation pressure at all and is shown for reference. The 
remaining lines correspond to the following cases: (a) solid - solar minimum, source region in the solar equator plane; (b) dotted 
- solar maximum, source region in the solar equatorial plane; (c) dash-dot - solar minimum, source region at ecliptic poles; (d) 
dashed - solar maximum, source region at ecliptic poles; (e) dash-dot-dot (practically overlapping with dash-dot) - force-free case, 
when straight-line, constant-speed motion is assumed. Right-hand panel: Differences in travel time between case (a) - (b) (dotted 
line); (a) - (c) (dash-dot); (b) - (d) (dashed); and (a) - (e) (dash-dot-dot). 



sensitivity band are up to 6 to 8 months, with equatorial atoms 
arriving more rapidly than the polar ones. 

The atoms that arrive simultaneously at 1 AU with energies 
within the IBEX sensitivity limit carry information from the 
source region that pertains to distinctly different epochs, span- 
ning a time interval larger than the solar-cycle period. 

3. Survival probabilities of H ENA 

3.1. Ionization model used 

The survival probability of H ENA against ionization Wion, be- 
tween the source region and the Earth orbit can be calculated 
from the formula: 



Wion = exp 



r 

•^'^l AU 



/3(t,r(t)) dt 



(4) 



where /3 (t, r (f)) is the ionization rate at time t along the trajec- 
tory r (f), ^source IS the time of "launch" of the atom in the source 
region, and f i au is the time of arrival of the atom at 1 AU; the 
integration goes backwards in time, which results in the absence 
of the minus sign in the exponent. 

The ionization rate fi is determined by all ionization pro- 
cesses affecting H ENA in the inner heliosphere, i.e. charge ex- 
change with solar wind particles (protons and alphas), collisional 
ionization by solar-wind electrons, an d ionization by solar p ho- 
tons, and was discussed extensively bv lBzowski et al.l (l2008l) for 
neutral interstellar gas in the inner heliosphere. The model of 
the ionization field relevant for heliospheric H ENA used in the 
present paper is discussed in greater detail in Appendix A; here, 
only the most relevant features are pointed out. 

To calculate the rate of charge exchange (Eq. IA.3b . we need 
the relative velocities between H ENA and solar wind protons. 
For trajectories considered in this paper, these velocities vary 
from an algebraic sum of a H ENA speed vena and the solar- 
wind speed \'sw at large distances from the Sun, Vrei - vsw + 
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Fig. 8. Adopted profiles of solar-wind speed (upper panel) and 
density at 1 AU (lower panel) for solar minimum and maximum 
conditions (respectively, solid and dotted lines). Corresponding 
profiles of the flux are shown in the lower panel of Fig. [17] 
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Vena, to a value appropriate for perihelion, v'rei 
The change in relative speed occurs over a period of a few 
months, just before reaching Earth orbit, when changes in the 
H ENA velocity vector also occur, as discussed in the preced- 
ing section. Hence, the rate of charge exchange during the most 
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Latitude profiles of electron ionization rate at 1 AU and 5 AU, 
solar minimum and maximum 
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Fig. 9. Latitudinal profiles of the normalized electron ionization 
rate R^/3ei{R) at solar minimum (solid lines) and maximum (dot- 
ted lines), calculated at 1 AU and at 5 AU. The 5 AU lines are 
the lower ones. 

Equatorial and polar electron ionization rate 
normalized to 1 AU 
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Fig. 10. Normalized radial profiles PeiiR) of equatorial (solid) 
and polar (dotted) rates of electron-impact ionization of the in- 
coming H ENA, defined, respectively, by Eq. (IA.13l l and ( IA.14l i. 



important time interval of perihelion approach, differs from the 
usually-adopted 

Since the solar wind is bi-modal, that is slow and dense in 
an equatorial band and fast and rarefied in the polar regions, 
the field of ionization of the incoming H ENA has latitudinal 
anisotropy, especially during solar minimum. The range of the 
slow-wind band expands from solar minimum to maximum with 
some phase shift between the north and south hemispheres, so 
that during solar maximum the slow wind persists at practically 
all latitudes. This fact is taken into account in the model used. 

Adopted profiles of solar-wind speed and density are defined 
by Eqs ( IA.4b . dA.Sb and ( IA.6b and shown in Fig. [8] The profile 
widths change in time according to Eq. (|A.6b . For solar mini- 
mum, the epoch 1995.0 was adopted, and for solar maximum, 
epoch 2001.0. The model features some north-south asymme- 
try, which agrees with the results of observations discussed in 
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Fig. 11. Solar wind velocity (upper panel) and density values 
(lower panel) from the OMNI-2 series, averaged by Carrington 
rotations. 
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Fig. 12. Carrington-averaged photoionization rate of H ENA 
during solar cycle based on SOLAR 2000. 



Appendix A.l. It was assumed that the solar-wind speed does 
not depend on distance from the Sun and that proton density de- 
clines as 

Since the physical state of electrons in the solar wind de- 
pends on the slow/fast regime, the electron-impact ionization 
rate features a characteristic latitudinal anisotropy, as shown in 
Fig. |9] The figure also illustrates the rapid decrease in the elec- 
tron ionization rate with heliocentric distance, due to cooling 
of solar-wind electrons. The reduction in the electron ionization 
rate with solar distance is more rapid than in the case of charge 
exchange (Fig.fTOli. The model of electron ionization rate, as pre- 
sented in greater detail in the Appendix, depends on the density 
of solar- wind protons at 1 AU nswiip, (Eq. ( |A.5b ). 
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Effect of radiation pressure on survival probability 
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Fig. 13. Survival probabilities of H ENA for a spherically- 
symmetric radiation pressure and ionization field. Solid line, 
left-hand vertical scale: w,„,„ for the solar Lyman-a flux equal 
to /„„■„, dotted line: w„a^^ for the solar flux equal to dash- 
dot, right-hand vertical scale: ratio of survival probabilities 



Both the density and speed of the solar wind show strong 
fluctuations (Fig. fTTT i. In the calculations, it was, however, as- 
sumed that the equatorial densities and velocities are constant 
(and their values can be found in Fig.|8]l; the effect of solar- wind 
fluctuations on survival probabilities is discussed separately. 

The photoionization rate also shows short-time fluctuations 
and distinct secular variations that are correlated with the solar- 
cycle phase (Fig.[T2]l. In the calculations, it was assumed that the 
photoionization rate is invariable and spherically symmetric; for 
solar minimum at 1 AU, [iph.mm - 0.8 x 10"^ s"\ and for solar 
maximum, ySpft mo r = 1.4 x 10"^ s"'. It also changes with solar 
distance, being proportional to 1 /r^. 

taruntman ( 1990) pointed out an additional, indirect channel 
of photoionization losses: a H atom is first excited from ground 
state to a higher state and then photoionized by solar photons 
that have a lower energy than needed to perform the ionization 
directly. This effect increases the photoionization rate by 10% at 
1 AU, but since its efficiency decreases with solar distance far 
more rapidly than it is reduced to just ~ 1% of the direct 
photoionization rate at 3 AU. Given the uncertainty in the pho- 
toionization rate, it can be neglected as an additional mechanism 
of photoionization losses in the incoming H ENA. 

3.2. Survival probability as function of radiation pressure 

The probability of survival of H ENA traveling to Earth orbit 
depends on a convolution of the local instantaneous ionization 
rate along their trajectories with the effects related to their kine- 
matics. Given a fixed ionization field, the probability of survival 
depends on the effective exposure time to the ionizing effects, 
which depends on the radiation pressure. This phenomenon is 
illustrated by the following example. 

When we assume that both the solar Lyman-ff flux and the 
ionization field are spherically symmetric and that the solar- wind 
parameters at all latitudes are similar to those at the solar equator 
(see Figs.[8]and|9]l, then the probability of survival of an H ENA, 



Role of vai'ious effects in net survival probability 
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Fig. 14. Contribution of various parameters to ionization losses 
of H ENA at Earth orbit. Shown are ratios Wi„„ Jwio„, where w,o„ 
is the probability of survival obtained with the use of the full 
ionization model including charge exchange ( with the cross sec- 
tion according to lLindsav & Stebbingsl ( l2005h ). electron ioniza- 
tion, and photoionization, calculated for trajectories in the solar 
equator plane, and vv',o„ .v is the probability computed along the 
same trajectories with photoionization switched off (diamonds), 
with electron ionization switched off (star s), and with a change 
of the cha rge exchange cross s ection from Lindsay & Stebbingsl 
( l2005h to lPhaneuf et alj (Il987l) (boxes). Solid fines correspond 
to the values relevant for solar minimum, broken lines to solar 
maximum. Additionally shown is the influence of the change of 
radiation pressure from solar minimum to maximum (triangles) 

^ ion,max / ^ ion,mm • 

reaching perihelion 1 AU, does not depend on the arrival direc- 
tion of the atom. It depends solely on its energy in the source re- 
gion and on the radiation pressure. This is illustrated in Fig.[T3] 
We show the survival probabilities of an H ENA which reaches 
perihelion at 1 AU, as a function of energy for solar Lyman-a 
fluxes /min and /max. They differ for energies Eiau < 0.05 keV, 
as demonstrated by the dash-dot line in the figure. For atoms 
of the lowest energies, the differences in survival probability re- 
sulting solely from to a change of /tot from /min to /max may be 
as high as ~ 40%. 

In reality, however, the ionization field is not spherically 
symmetric, and another imprint of the atom dynamics might be 
expected due to the latitudinal anisotropy of the solar Lyman- 
a flux. A modulation of the H ENA survival probability as a 
function of latitude of the approach direction would be expected. 
Calculations (not shown) suggest, however, that this effect yields 
a modulation with an amplitude equal to just 5% of the average 
value for the lowest energies accessible by IBEX; this amplitude 
drops down below 1% for E\a\] = 0.03 keV, so this effect is 
negligible. 

3.3. Contributions of various loss channels 

Contributions from various ionization channels to the net sur- 
vival probability of H ENA at Earth orbit are presented in 
Fig. [14] Shown are ratios vv,on,A/w,o„ calculated for various tra- 
jectories, where w,o„ is the survival probability, obtained from a 
full model taking into account all three ionization channels, and 
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Probabilities of H ENA survival normalized to ecliptic 





Fig. 15. Latitudinal profiles of ratios of H ENA survival probabilities w,o„ {(p) Iwum (0) for various H ENA etiergies at Earth orbit, 
shown for solar minimum conditions and for 4 locations at Earth orbit. The energies are indicated in the figure; the profiles without 
energy labels are for (from top to bottom) 0.208 keV, 0.403 keV, 0.781 keV, and 5.9 keV. The absolute values of probabilities for the 
reference atoms, arriving from the ecliptic plane, are indicated in Fig. [13] Shown are results of 360°scans in the planes perpendicular 
to the local antisolar vector, parametrized by the ecliptic latitude of the aiTival directions of the atoms. 



Wion,x is the probability of survival obtained when a factor "x" 
is switched off. The probabilities were calculated for trajectories 
contained in the solar equatorial plane, i.e. fully immersed in the 
slow solar wind. In such circumstances, the most intensive ioniz- 
ing factor after charge exchange with solar-wind protons is pho- 
toionization, whose share is a strong function of atom energy. 
In the case of the slowest atoms (~ 0.01 keV), photioinization 
contributes as much as ~ 50% of the losses, but only ~ 5% for 
atoms with isiAu =6 keV. During solar minimum, when polar 
regions are engulfed by the fast wind and the charge exchange 
rate is much lower than in the slow wind, photoionization how- 
ever, can be a dominant source of losses of H ENA originating 
in polar latitudes. 

Electron-impact ionization produces up to ~ 20% of the ion- 
ization losses of slowest H ENA in the slow solar wind and 
its contribution decreases to below 10% at /iiAu - 0.05 keV. 
During solar minimum the electron ionization rate is far more 
anisotropic than the charge exchange rate (cf Fig.|9]and the lower 
panel of Fig. [TTI i; we therefore expect that the contribution of 



electron ionization to the losses of H ENA approaching Earth 
orbit from polar latitudes during solar minimum will be lower 
than in the case of atoms originating in the equatorial latitudes. 

The contributions of photoionization and electron ionization 
to the net losses of H ENA depend weakly on radiation pressure 
(the change from solar minimum to solar maximum is just a few 
percent), as illustrated by the dotted lines in Fig. [14] 

In the case of the slowest H ENA at Earth orbit, the most 
significant influence on the change in survival probability be- 
tween solar minimum and maximum is the change in radiation 
pressure. This is illustrated by the line of triangles in Fig. [14] 
The amplitude of change in survival probability of H ENA at 
transition from solar maximum to solar minimum decreases 
to below the level of the contribution of photoionization for 
£^iAU - 0.015 keV, and below the contribution of electron ion- 
ization for E\ Au - 0.03 keV. 

The uncertainty in the cross section for charge exchange 
(boxes in Fig. [I4] | produces an uncertainty of ~ 10% in the es- 
timate of the survival probability of H ENA at the lower end of 
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south polar arrival direction to minimum value 
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Fig. 16. Upper panel: Ratios of survival probabilities of H ENA 
approaching Earth orbit from the south ecliptic pole to the low- 
est probability along the Earth orbit for various H ENA energies 
at 1 AU: from top to bottom 0.010 keV (solid line), 0.015 keV 
(dots), 0.029 keV (dash-dot), 0.107 keV (dashed), 0.209 keV 
(dash-dot-dot), and 5.9 keV (solid). Lower panel: Ratios of the 
highest to lowest probability of survival of H ENA originating in 
the south ecliptic pole region as a function of energy at Earth or- 
bit for solar minimum (solid) and maximum conditions (dotted 
line). 



energies at Earth orbit, and decreases, with an increase of energy, 
to the level of a few percent at ~ 6 keV. 

3.4. Anisotropy of survival probability 

The probability of survival of H ENA reaching Earth orbit at 
perihelia of their trajectories exhibits an anisotropy as a function 
of the latitude of arrival direction, which is a result of combi- 
nation of the anisotropy of the ionization rate and of geometry 
effects due to the non-vanishing inclination of Earth orbit to so- 
lar equator. In the adopted model, the ionization field has an ap- 
proximate planar symmetry with respect to the solar equator and 
IBEX will carry out observations remaining bound in the eclip- 
tic plane, which is inclined to the equator at an angle ~ 7.25°. 
The anisotropy is illustrated in Fig.fTSl 

Assuming the anisotropic ionization rate as shown in Figs. [8] 
and|9]and an invariable, anisotropic solar Lyman-o' flux (Eq. (O), 
the distributions of the survival probability of H ENA of a given 
energy, reaching Earth orbit at ecliptic longitude 164° (which is 
the crosswind point and simultaneously the point in the eclip- 
tic with lowest heliographic latitude -7.25°; lower-left panel) 
depend strongly on the latitude of aiTival direction. The pole-to- 
ecliptic probability ratio for the lowest-energetic atoms exceeds 
2, but decreases with energy increase, reaching ~ 1.1 for the 
highest-energy atoms. Latitudinal profiles of the survival prob- 
ability show some north-south asymmetry related to the north- 
south asymmetry present in the ionization model, but there is no 
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Fig. 17. Upper panel: Probabilities of survival of H ENA ap- 
proaching Earth orbit at ecliptic longitude 344° as a function of 
ecliptic latitude of the arrival direction for 1 AU energies equal to 
0. 107 keV (the lower group of lines) and to 0.209 keV (the upper 
group of lines) for three phases of solar cycle (from top to bottom 
in both groups of lines): minimum, intermediate and maximum 
(respectively: solid, dots, dash). Lower panel: Adopted solar- 
wind fluxes as a function of heliolatitude for the three phases 
of solar cycle (minimum: solid, intermediate: dots, maximum: 
dash), used to calculate the survival probabilities of the H ENA 
at Earth orbit shown in the upper panel. The geometry of the up- 
per panel corresponds to the geometry of the lower-right panel 
of Fig. [B] 



asymmetry in ecliptic longitude (there exists a left-right symme- 
try). 

A change in the observation point along Earth orbit produces 
a distortion of this picture, which is entirely due to the shift in 
the observation point with respect to the solar equator. This effect 
is shown in the upper-left panel of Fig. [15] (a point in the solar 
equator plane and simultaneously the upwind point of Earth or- 
bit), where one can see - for the ionization field unchanged - an 
appreciable increase in the north-south asymmetry as compared 
with the previous case and appearance of an asymmetry in eclip- 
tic longitude (the lines in the figure are hysteresis-like, there is no 
left-right symmetry). Another change in the observation point to 
the other crosswind point (Acci = 344°, a location with the high- 
est heliolatitude, equal to -1-7.25°) produces another change in the 
relative survival probability profiles (lower-right panel): because 
of the shift with respect to the solar equator, the north-south 
asymmetry is strengthened considerably and the left-right asym- 
metry disappears. In the case of the slowest atoms, the north 
pole-to-ecliptic survival probability ratio appreciably exceeds 3. 
A transition to the downwind point of the Earth orbit (the other 
node of solar equator; upper-right panel) of course restores the 
picture from the upwind point, but the local left and right direc- 
tions are swapped. 

Changes in the survival probabilities of H ENA originating 
in polar regions are shown in Fig. [161 The upper panel of the 
figure presents the ratio of the survival probabilities wion^s i^eci) 
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of H ENA originating in the south ecliptic pole and reaching 
Earth orbit at ecliptic longitude Aeci, to the lowest probability 
along the Earth orbit Wio„,s (min). The ratios are shown for var- 
ious energies at Earth orbit during solar activity minimum. The 
changes in survival probability shown are solely due to geomet- 
ric effects. For atoms of the lowest energies at Earth orbit, the 
yearly amplitude of the ratio w,o„ 5 (Aed) / w,on,s (min) is approx- 
imately 30% and decreases rapidly enough with increasing en- 
ergy. For energies that exceed ~ 0. 1 keV, it decreases below 10% 
(Fig. [16] lower panel). The maximum of survival probability of 
the H ENA originating in the south ecliptic pole and observed by 
IBEX will occur at ecliptic longitude ~ 164°, i.e. at the transition 
from February to March, and in the case of atoms originating in 
the north ecliptic pole, the course of survival probability as a 
function of ecliptic longitude will have a similar shape, but with 
a phase shift of 180° (half a year). 

It appears that anisotropies larger than those presented 
should not be expected. The evolution in the anisotropy of H 
ENA survival probability for three phases of solar cycle at A^ci - 
344° for £1 Au - 0.1 keV and ~ 0.2 keV is shown in Fig. [T71 
along with the corresponding profiles of the solar-wind flux. In 
addition to the solar minimum and maximum phases, as previ- 
ously, a phase of intermediate solar activity was added, when the 
band of slow solar wind had expanded to mid-latitudes. For this 
phase the solar Lyman-a flux /tot = 4.7 x 10" cm"^ s"' and the 
photoionization rate /3pi, = 1.1 x 10~^ s"' were adopted. The lat- 
itudinal profiles of the solar wind flux for these phases of solar 
activity are shown in the lower panel of Fig. [T7] The expansion 
of the latitudinal range of slow solar wind, related to the increase 
of solar activity, produces a flattening and widening of the pro- 
files of Wion (0) / Wio„ (0), with an almost complete disappearance 
of anisotropy at solar maximum. 

3.5. Fluctuations of survival probability 

Because of strong fluctuations in both the solar wind (Fig. fTTT i 
and the photoionization rate (Fig. [T2b . one should expect ap- 
preciable fluctuations in the survival probabilities of H ENA at 
Earth orbit. Available data enable us to estimate these fluctua- 
tions for arrival directions contained in the ecliptic plane, but 
because of insufficient data the estimate for higher latitudes will 
be less accurate. Since Ulysses observations imply that the fast 
wind, which is characteristic of high heliolatitudes, is aff'ected 
by smaller fluctuations than the slow wind, we expect that fluc- 
tuations for H ENA originating in the polar regions will be lower 
than for the ecliptic directions. 

Based on a simplified kinematic model (with a force-free 
motion of H ENA), the ratios of the H ENA survival probabili- 
ties at Earth orbit w,+i /w,- were studied for various energies. For 
solar-wind conditions that are averaged by i-th Carrington ro- 
tation, w, is the probability of survival. Solar-wind parameters 
from the OMNI-2 time series were used (Fig. fTTT i. An example 
histogram of these ratios for Eiau = 0.1 keV is shown in the 
upper panel of Fig. [T8j the distribution is close to normal. The 
dependence of the standard deviations of these ratios / w, on 
energy, CTstdoev (wi+\/wi, E\ au), is illustrated by the dotted sym- 
bols in the lower panel of Fig.[T8j This dependence can be ap- 
proximated by the formula: 

o-stdDev (w/+i /w/, £1 Au) = 0.075 cxp [-0.03 In^ (£1 au)] ^Tau ' 

(5) 

with energy Ex au expressed in keV, and is shown in the lower 
panel of the figure with the solid line. This suggests that the 
fluctuations 5w m the survival probability at Earth orbit, due to 
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Fig. 18. Upper panel: An example histogram of fluctuations of 
quantity w,+i/vvi, where w, is the probability of survival of H 
ENA calculated with solar-wind equatorial parameters averaged 
by Carrington rotation periods for an H ENA energy at Earth 
orbit equal to 0.1 keV. The solar- wind parameters were taken 
as the averages of appropriate intervals of the OMNI-2 time se- 
ries, shown in Fig. [TT] Lower panel: Standard deviation of the 
monthly survival probabilities ratio w,+i/w, as a function of H 
ENA energy at Earth orbit (dots) and a plot of the approxima- 
tion formula specified in Eq. (|5]l. 

fluctuations in the ionizing factors, should range 0.2 Wexp, where 
We.vp is the expected probability value at the lower end of the 
IBEX energy band, to negligible values {6w - Q'm the case of 
H ENA from the high limit of this band. For H ENA of energies 
expected from the heliosheath 0.1- 0.2 ke V, the fluctuations are 
on the order of ~ 0.2 Wexp - 0.15 Wexp, which is smaller than the 
amplitude of the latutidinal modulation related to the solar-wind 
anisotropy (Figs.[T5]and[T7]l. 

4. Summary and conclusions 

We investigate modifications of H ENA energy at Earth orbit 
with respect to their energy at 100 AU from the Sun and their 
probabilities of survival against losses induced by local helio- 
spheric ionization processes. The range of energy at Earth or- 
bit that we have studied corresponds to the sensitivity range of 
the planned NASA SMEX mission IBEX, with emphasis on the 
H ENA energies expected from the inner heliosheath (mostly 
0.1 - 0.2 keV, up to ~ 1 keV). Realistic, observation-based 
models of both ionization processes and radiation pressure have 
been used. The ionization model reproduces the time-dependent 
spatial anisotropy of the solar wind, and the radiation-pressure 
model takes into account the dependence of ju on both the radial 
velocities of the atoms and on the heliolatitude (Figs.fTllSllQlfTTl 
and[T2and Eqs dD, ©, and Q). 

Based on the results of our numerical simulations, we pro- 
pose that, due to radiation pressure, H ENA reach Earth orbit 
without a significant changes in energy and direction apart from 
the atoms with energies below 0. 1 keV during high solar activ- 
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ity (Figs, m |6]l. The probability of H ENA survival is a function 
of the solar-cycle phase, the ecliptic longitude of the location at 
Earth orbit, and the ecliptic latitude of the approach direction. 
Especially important in this respect is the latitudinal range of 
the slow solar-wind region (Figs. [TSjandfTTli. For all cases con- 
sidered, the time of flight from the source region at 100 AU to 
Earth orbit is equal to the force-free travel time, to an accuracy 
of 2 solar rotation periods, down to 1 AU energies ~ 0.05 keV. 

The survival probability of H ENA declines rapidly with de- 
creasing energy at Earth orbit from ~ 0.8 for £1 au = 6 keV to 
~ 0.4 for £1 AU = 0.2 keV and ~ 0.3 for £iau = 0.1 keV for H 
ENA approach directions within the ecliptic plane (Fig. [13). A 
yearly modulation of the survival probability of H ENA originat- 
ing in the polar regions, is expected due to the yearly movement 
of the Earth-bound detector in heliolatitude. The yearly ampli- 
tude of the maximum-to-minimum probability ratio decreases 
with increasing energy at Earth orbit; for Eimj - 0.1 keV the 
ratio is equal to ~ 1.15, and for 0.2 keV to ~ 1.1 (Fig.fTSIl. The 
courses of the modulation for the north and south poles are in 
antiphase. 

Apart from a short time interval at solar maximum, one 
also expects an appreciable modulation of the survival proba- 
bility as a function of heliolatitude (Figs. [15] and [TTT i. During 
solar minimum in the case of H ENA with E[au = 0.1 keV 
the amplitude of this modulation should not exceed 1.6, 1.3 for 
EiAU = 0.2 keV and 1.1 for /iiAu = 6 keV. Appreciable left- 
right and north-south asymmetries are expected too, changing 
with position at Earth orbit (Figs.[T5]and[T7li. as well as a distinct 
correlation of the survival probability profiles with the latitude 
structure of solar wind. The amplitude of the latter modulation 
decreases with widening of the slow solar wind region (Fig.fTTli. 

Because of fluctuations of the solar wind and of solar EUV 
radiation at time scales shorter than solar rotation period fluctu- 
ations of the H ENA survival probability at Earth orbit are ex- 
pected. Based on the solar wind observations at 1 AU and on 
the simulations discussed in the paper a phenomenology func- 
tion expressing the amplitude of these fluctuations with respect 
to the survival probability as a function of the H ENA energy 
at Earth has been found (Eq.©): it decreases with an increase 
of energy and becomes lower than 15% for Eimj - 0.1 keV 

(Fig.m- 

We have also studied the contributions of various factors af- 
fecting the losses of H ENA traveling from the source region 
to Earth orbit. The uncertainty in the charge-exchange cross- 
section yields uncertainty in the survival probability approxi- 
mately of a few percent, much lower than the amplitude of the 
probability fluctuations due to fluctuations in the solar wind. 
The strength of radiation pressure affects the survival probabil- 
ity on the level of ~ 20% for £1 au - 0.1 keV and -15% for 
EiAu - 0.2 keV. Photoionization contributes 20% and 15% of 
the losses at the energies /siau - 0.1 keV and £1 au - 0.2 keV, 
respectively, and far more for approach directions from high lati- 
tudes. Electron ionization for energies above 0. 1 keV contributes 
less than 10% of the losses (Fig.[T4li. 

All modulation factors discussed increase appreciably - al- 
though to varying degrees - with a decrease in H ENA energy at 
Earth orbit. They should, however, disappear during solar max- 
imum, apart from those related to solar-wind fluctuations, when 
the large-scale structure of the solar wind is spherically symmet- 
ric. 

The analysis presented in this paper has neglected the obser- 
vational effect of a proper motion of the detector, which will be 
the subject of future studies. Based on results obtained thus far it 
is concluded that the most important elements needed for a suc- 
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Fig.A.l. Cross sections for charge exchange between H ENA 
and protons (solid line) and alpha particles (dotted line) for col- 
lision speeds typical of the inner heliosphere. Broken line is the 
alpha charge exchange rate multiplied by a typical abundance of 
alphas with respect to protons, which is assumed to be 0.04. 

cessful interpretation of H ENA measurements, such as those 
planned for the forthcoming IBEX mission, are monitoring of 
the local solar wind, the latutudinal structure of the solar wind, 
and the solar Lyman-or flux. The solar EUV flux in the spectral 
range responsible for photoionization of hydrogen is also im- 
portant, and especially its dependence on heliolatitude. Taking 
into account the local effects of H ENA transport from source 
to Earth orbit, we expect that IBEX should be able to discover 
departures from symmetry in the heliospheric interface of ap- 
proximately a few percent. 

Appendix A: Ionization processes destroying 
heliospheric H ENA 

A.1. Charge exchange with solar-wind protons and alpha 
particles 

The local rate of charge exchange between a H ENA traveling 
with velocity Vena and a proton or alpha particle traveling with 
velocity Vp and a member of solar-wind population described by 

a local distribution function fp (vp^ is defined by the formula: 

Pcx = J o"cx (|v£iv^ - Vp|) \vena - Vp\ fp (vp) dVp, (A. 1) 

where cr^x is the charge-exchange cross-section, which is a func- 
tion of the relative velocity between the colliding particles Viei = 
Vena - Vp- Since the solar wind is hypersonic, the thermal spread 
of its positively-charged particles can be neglected. The distribu- 
tion function fp is then replaced with a delta function and with 
the substitution 

Vrel = IVrell = \Vh - VswI , (A.2) 

where vsw is the local solar-wind velocity vector, Eq. dA.U sim- 
plifies to the form: 

yScx = O-cx (Vrel) v.^i «sw (A.3) 
where «sw is the local density of solar wind protons. 
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Since the solar-wind velocity in the inner heliosphere can 
be treated as purely radial and the incoming H ENA travel 
almost radially towards the Sun, we can determine approxi- 
mate limits of the relative velocities between the colliding solar- 
wind ions and the H ENA that at 1 AU have energies within 
the IBEX sensitivity limits: from ~ 400 km/s for the slow- 
est ENA and slow solar wind to ~ 2000 km/s for the fastest 
ENA and fast solar wind. This allows us to assess the cross sec- 
tions for the charge exchange of H ENA with solar-wind pro- 
tons and alpha particles: for a H H- He^^ collision the cross sec- 
tion varies from 1.3 x 10"'* cm^ for the low-energy collision 
to 10 cm^ for the high-energy collision (iPhaneuf et al.lll987l) 
and for the H + collision respectively fro m 1.8 x 10"'^ cm^ 
to 5 X 10"= cm^ dLindsav & Stebbingsl2005l) . The cross sections 
for charge exchange with protons are larger than the cross sec- 
tions for charge exchange with alphas only for relative velocities 
below ~ 1000 km/s (Fig. lA.lb : since the abundance of solar- 
wind alpha particles wi th respect to protons is , however, typi- 
cally just a few percent dvon Steiger et al.ll2000l) . the net contri- 
bution of the charge exchange with alphas is negligible for slow 
H ENA and up to 10% for the fast wind and fastest ENA. 

The precise value of the charge exchange cross-section 
has been a subject of debate. Fortunate ly for the topic of 
this paper, estimates by v arious authors dBarnett et al.l 119901: 
iLindsav & Stebbingsl l2005l) appear to have converged within 
~ 10% for the collision speeds relevant to H ENA charge ex- 
change with solar wind protons (apart from the highest energy 
wing, for which they increa se to ~ 20%). M ore details about 
this issue can be found in iFahr et akl d2007l) . In the simula- 
tions discussed in this paper, the cross section recommended by 
iLindsav & Steb bings ( 2005) was used. 

The solar wind is known to be anisotropic, consisting of an 
equatorial band of slower expansion of a denser plasma (inter- 
leaved with high-speed streams) and polar caps of a steady, rar- 
efied, but much faster wind. The latitudinal range of the jagged 
boundaries of the slow wind band vary during the solar cycle, 
being closest to the equator during solar minimum and expand- 
ing gradually towards the poles with an increasing level of solar 
activity, before engulfing the poles during the solar-maximum 
phase. The north and south hemispheres appear to be shifted in 
phase by a year or two. 

Details of this evolution are unknown. Latitudinal structure 
of the solar- wind was measured in situ by Ulysses, mostly during 
fast l a titude scans (Phillips et al. 1995; McComas et al. 1999, 
I2OOOI I2OO3I) . These measurements, although highly accurate, 
provide only point-like data on density and velocity without a 
global snapshot picture. In particular, they provide little informa- 
tion about the longitudinal pattern of the fast/slow wind bound- 
ary. 

Another technique to study the solar wind structure is map- 
ping the solar wind indirectly, by analysis of the solar Lyman- 
a radiation backsca ttered off the heliospheric neutral hydrogen 
dBertaux et al.lll99"5h . The dense and slow equatorial solar wind 
is able to ionize the neutral interstellar hydrogen inside the helio- 
sphere more efficiently than the fast but rarefied wind in the polar 
regions, and thus carves a trough in the gas distribution, seen as a 
dimmer band in the global map s of heliosp heric glow, referred to 
as the heliospher i c groo ve (e.g. lKvrola et al.lll998l:lBertaux et al.l 
Il999l) . iBzowskil (|2003|) proposed a quantitative model to infer 
the latitudinal span of the equatorial band of the slow wind from 
the observations of the latitudinal range of the groove, employed 
by Bzowski et al. (2003) to determine the evolution of the slow- 
wind region from solar min imum in 1996 to solar maximum in 
2001 ■ iBzowski et aTl (l2008l) pointed out that there is a correla- 



tion between the latitudinal span of coro nal holes and the bound- 
aries of the slow wind region. Similarlv. 'Ouemerais et al.l(l2007h 
demonstrated a correlation between coronal features observed in 
white light and the structure of the solar wind. Based on such 
correla tions and on remote-sensing information, iBzowski et al.l 
(I2008h suggested a phenomenology, proxy-based model for the 
evolution of structure in the solar wind, which forms the basis 
of construct a of solar-wind model constructed to calculate the 
charge-exchange losses discussed in this paper 

The Lyman-Q- remote-sensing technique is only able to re- 
cover the structure in the ionization rate of the low-energy helio- 
spheric hydrogen, which is correlated with the solar-wind flux. 
Thus, it is able to yield a global structure of the solar wind over 
an extended interval of time, but is unable to recover separately 
the solar-wind speed and velocity, which are essential to cal- 
culate the charge exchange rate of the incoming H ENA (see 
Eq. dAJt and dlS). 

In addition to Ulysses in situ and SWAN remote-sensing ob- 
servations we, fortunately, also have continuous monitoring of 
the solar wind in the ecUptic plane. The data retrieved from var- 
ious experiments were cross-calibrated and c ompiled into the 
OMNI time series (iKing & Papitashvliill2005h . which covers a 
few solar cycles until present (Fig. fTTT l. As one can see in the 
figure, effects of the level of solar activity are not so obvious as 
they were in the case of the solar EUV output shown in Figs.[T] 
and [12] We note a downward trend in solar- wind density, which 
initiated about the previous solar minimum and appears to con- 
tinue until the present; a clear, periodic signal appears, however, 
to be missing. Velocity shows some correlation with solar ac- 
tivity, increasing during solar maximum and decreasing during 
solar minimum. The amplitude of these variations, however, is 
only about 10%, so the flux is mostly governed by the density 
changes, and decreases steadily during the past solar cycle. 

With these data sets and observations in hand, we attempt to 
construct an evolutionary model of latitudinal density and veloc- 
ity structure of the solar wind. We develop the following model 
by assuming that (i) the values of the fast solar-wind density 
and velocity at the poles are invariable (although the north- and 
south-pole quantities may differ from each other), that (ii) the 
fast/slow wind boundaries are parallel to the solar equator (not 
exactly true, but close to the truth when averaged over all solar 
longitudes) and equal to the values measured by Ulysses during 
the fast latitude scan, and that (iii) the equatorial values corre- 
spond to the values provided by the OMNI series (which agree 
with Ulysses values at its passages through ecliptic), and adopt- 
ing the variations in the north and south boundaries of the slow 
wind band as inferred from the SWAN Lyman-a mapping and 
the coronal proxies. 

The latitudinal profile of the solar-wind speed is defined by 
the formula: 

I'SW.pol) 



(vsw.pol + 0) + (vsw.eqtr (0 

t>N it) - 4>s (0 



(A.4) 



exp 



-ln2 



(pN(t)-^sit) 



where 4> is heliographic latitude and is a shape factor, which in 
the present study was adopted to be = 8; vsw.poi is the average 
solar-wind velocity at the poles and the term (vsw.poi + 4>) de- 
scribes the north-south asymmetry of the polar velocities; the 
term (vsw.poi + 5y </>) + (vsw.eqti- (0 - vsw.poi) for = corre- 
sponds to the solar wind speed at solar equator; and the term 

- In 2 (j^^^^z^^) describes the latitudinal dependence of 
the velocity. 



exp 
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The density profile nsw (tp, is defined by a similar formula: 

nsw (<!>, t) - (nsw.pol + Sn (/>) + ("SW.eqtr (0 - "SW.pol) (A.5) 



X exp 



-ln2 



' (0 - <f>s (0 



( t) - 4>S i t) 



where R - r/rE (fE = 1 AU). The parameters in Eqs 
(IA.4b and (IA.5l l were fitted to the Ulysses data from the first 
fast latitude scan and the results are the following: vsw.poi - 
766.0 km/s, = 0.14822 km/s/deg, nsw.poi = 2.5180 cm^^ 



2.0009 ■ 10 cm ^/deg. The evolution of the north and 



sout h boundaries of the eq uatorial slow-wind region is adopted 
from lBzowski et all (l2008h : 

(pN,S (0 = <f>ON,S + <f>l N/S exp [- COS^ (w0iv_s ?)] , (A.6) 

where ^o^ai = 36.4°, - -22.2°, lo^^n - 0.58251, 0o,s - 
-40.600°, = 20.0°, = 0.58226. ForreaHstic modeling, 
the values of parameters neqtr, and Veqtr for a given time t (which 
must be expressed in decimal years) should be adopted from the 
OMNI time series (Fig. [TTI ): in the fits to the first Fast Latitude 
Scan by Ulysses used throughout this paper, they were equal to 
7.02 cm"^ and 426.7 km/s, respectively. The density and velocity 
profiles returned by the model for solar minimum and maximum 
are shown in Fig. [8] and the flux in the lower panel of Fig. [17] 



A.2. Photoionization 

The rate of photoionization of heUospheric hydrogen is surpris- 
ingly poorly understood. To measure the value /3ph, we must in- 
tegrate solar spectral flux F (A) in the energy range above the H 
photoionization threshold at 13.59 eV (Aq = 91.175 nm), mul- 
tiplied with the energy-specific photoionization cross-section 
o-ph (A): 



ySph^ r F(A) <T^i,(A)dA 



(A.7) 



Details can be found, e.g., in 'Qga wa et al.l (Il995h . The solar 
spectrum in the relevant frequency range consists of a contin- 
uum, which decreases rapidly with decreasing wavelength, and 
of a flickering line-component, which is responsible mostly for 
solar cycle-related variations in the net photoionization rate. 
Hence, monitoring the neutral H photoionization rate requires 
a continuous monitoring of the solar spectrum shortwards of the 
photoionization threshold. While such measurements were per- 
formed by SOHO, there exists a spectral gap in the coverage 
just short of the ionization threshold, which must be filled with 
proxies. In the past, a very simple proxy scheme based on the so- 
lar 10 .7 cm radio flux was in use (iRuciiiski et al.,i.l996: Bzowskil 
I2001h . Recently, a more ingenuous proxy scheme \ yas developed 
within the framework of the SOLAR 2000 model (iTobiska et all 
[2000), which was adopted in the present paper 

The photoionization rate shows a strong correlation with the 
solar Lyman-a flux. It varies distinctly during the solar cycle, 
with a mean value equal to about 1.1 x 10"^ s"' and an ampli- 
tude of approximately 30% (Fig.[T2li. Thus, the photoionization 
rate is equal to about 20% of the equatorial charge exchange 
rate. It ha d been often a ssumed to be spherically symmetric, but 
studies bv lWittd (12004 . concerning the photoionization rate o f 
neutral heliospheric helium, supported bv lAuchere et al.l (l2005l) . 
imply that the hydrogen photoionization rate should HAVE a he- 
liolatitude anisotropy, which varies during the solar cycle. Due 
to insuflicient data, this efifect has been omitted from our present 
analysis. 



Since the inner heliosphere is almost perfectly transparent to 
solar EUV radiation of wavelengths shorter than the hydrogen 
ionization threshold, the H photoionization rate decreases with 
solar distance as l/r^. 



A.3. Electron impact ionization 

Construction of a model of electron-impact ionization in the in- 
coming H ENA is hampered by the fact that the behavior of 
the electron fluid in the solar wind is not fully understood. It 
is known that their local distribution function can be decom- 
posed into 3 components: a warm core, a hot halo (both approx- 
imated by Maxwellians), and a fluctuating spur, stretched along 
the local magnetic-field direction (Pilipp et al. 1987). Electron 
density can be taken from quasi-neutrality and continuity con- 
ditions of the solar wind and hence can be adopted as equal to 
the local proton density -H 2 x alpha density. The temperature 
in the region of interest (i.e. outside 1 AU at all heliolatitudes) 
is more uncertain. From measurements performed by Ulysses, 
it appears that it differs from that of protons and that its radial 
profile can be described by power laws with exponents in the 
range between isothermal and adiabatic values. The profiles in 
the fast wind differ from those in the slow wind. The content of 
the halo population is a function of heliocentric distance and of 
the sola r- wind regime (fast / slow). These findings were recapitu- 
lated by* Maksimovic etaP (l2000l) . The measurement of the fol- 
lowing quantities are, however, unreliable: the exponents of the 
temperature power laws, the balance between the core and halo 
populations, and their evolution during the solar cycle. The con- 
struction of an accurate model of electron ionization is therefore 
presently impossible. Fortunately, electron ionization is a minor 
loss process in the incoming H ENA and thus an approximate 
model is sufficient. 

The electron ionization rate of an incoming ENA can be cal- 
culated from the formula: 



>6e 



f 



fe (V'rel) 0"^ (Viel) Vrel d Vrel 



(A.8) 



where CT^is the cross section for electron-impact ionization 
(lLotzlll967h . f, is the local electron distribution function, Vrei 
is the relative speed between the incoming ENA and a solar- 
wind electron, and v'lin, is the velocity corresponding to the ki- 
netic energy which is equal to the hydrogen ionization potential 
13.59 eV. At electron temperatures of the order of 10"^ - 10^ K, 
the thermal spread in the electron distribution function is large 
in comparison with the relative velocities between the incom- 
ing H ENAs and individual electrons, and v,f/ in Eq. (IA.8I ) can 
be superseded by the electron specific velocity v. Following 
iRuciriski & Fahj (il989l) . the electron-impact ionization rate can 
be calculated by assuming that the local distribution function of 
electrons is composed of two Maxwellian functions (core and 
halo), parametrized respectively by core and halo densities tic 
and nil and core and halo temperatures Tc and T/,. 

These calculations were performed independently for the 
fast- and slow -wind regimes usin g parameters for the slow wind 
adopted from lScime et al.l (11994 : 



Tc = 1.3 • lO-V"" **^ 
Th = 9.2 • 10^ r""-^*^ 
^ch = «/,/«.. = 0.06 r-"-25 



(A.9) 
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and, for the fast-wind re gime, from llssautier et alj (Il998l) and 
iMaksimovic et all (l2000h : 

T, = 7.5- lO^r-O'^ 
ThIT, = 13.57 (A. 10) 

U = 0.03 

In both solar wind regimes, the core and halo densities tic and 
are calculated from the equations: 

1 +2^„ 
tic = rip 
1 + gc/i 

rih = ^chric (A. 11) 

with the alpha abundance ^„ = 0.04, which is identical in both 
fast and slow wind regimes, ^ch is the halo-to-core density ra- 
tio defined for slow and fast wind in Eq. ( |A.9t and ( lA.lOl ). re- 
spectively, and Hp is the proton density adopted from the density 
model discussed earlier in this paper. 

Based on these equations, calculations of the electron ioniza- 
tion rate as a function of heliocentric distance were performed 
(separately for the fast and slow modes) and the results were ap- 
proximated by the formula: 

, / X [ c (XnRf- + b\\\R + a 

R Bei[R,nr,] - Hp e\p 5 r . 

\ [gilnRf + filnRf +elnR + d 

(A. 12) 

The corresponding radial profiles of the electron-impact rate for 
the slow and fast wind are presented in a numerically-optimized 
form in Eq. ( IA.13I ) and ( IA.14I ). respectively: 

ln/?(541.691n/?- 1061.32)+ 1584.32 
!!p.^(\nR- 29.11) ((\nR - 2.02) in + 2.^(^^,13) 



«2 



In/; (348.73 \nR - 917.39) + 2138.05 



P.i,f{R,np) = ^e0n^-18.97)((ln/;-2.53)ln/; + 5.7(^.14) 

As is evident in these formulae, the electron-impact ionization 
rates are parametrized by local proton densities Hp normalized 
to 1 AU and are fixed functions of heliocentric distance, which 
differ appreciably from the 1 /R^ profiles that are typical of the 
solar wind flux and photoionization rate, as shown in Fig.fTOl 

The global time- and latitude-dependent model of the elec- 
tron impact ionization rate is constructed by analogy with the 
proton-density model: 



el,eqtr -/^el.pol) 
- (pN it) - (pS it) 



(A. 15) 



exp 



-ln2 



(f>N it) - (f>S it) 



where (p^ it) , (ps it) are defined in Eq.( |A.6l ). and 

jSel,eqti- = ySel.s {R, Wcqti- (O) 

j6ei,poi = y8ei,f(/?,npoi(f)) (A. 16) 

= /3,l,tiR,6n) 

with solar-wind related densities neqti, W poi ta ken from the model 
of solar-wind density described by Eq. ( IA.5I ). and the models of 
the slow and fast wind, indicated by Eq. ( IA.13I ) and Eq.( IA.14l ). 
respectively, evaluated for Hp = ngqtr and Hp = Wpoi. 

Example latitudinal profiles of the electron ionization rate, 
normalized to 1 AU by R^ and calculated for 1 AU and 5 AU, 
are shown in Fig.|9]for solar minimum and maximum conditions. 
The figure illustrates the rapid decrease in the electron ionization 
rate with increasing heliocentric distance and the evolution in 
latitudinal shape of the rate during the solar cycle. 



Acknowledgements. This study is partly based on results of the International 
Space Science Institute Team on "Effects of Heliospheric Breathing Due to Solar 
Cycle Variations from Back-Scattered Ly-Alpha", whose support and friend- 
liness is gratefully acknowledged. The SOLAR2000 Research Grade histori- 
cal irradiances are provided courtesy of W. Kent Tobiska and SpaceWx.com. 
These historical irradiances have been developed with funding from the NASA 
UARS, TIMED, and SOHO missions. The OMNI data were obtained from 
the GSFC/SPDF OMNIWeb interface at http://omniweb.gsfc.nasa.gov This re- 
search was supported by the PoUsh MSRiT grants 1P03D00927 i N522 002 
31/0902. 



References 

Auchere, F. 2005, Apl, 622, 737 

Auchere, F, McMullin, D. R., Cook, 1. W., & et al. 2005, in ESA SP-592: Solar 
Wind 11/SOHO 16, Connecting Sun and Heliosphere, 327-329 

Bamett, C. F, Hunter, H. T., Kirkpatrick, M. I., et al. 1990, Atomic data for 
fusion. Collisions of H, H2, He and Li atoms and ions with atoms and 
molecules. Vol. ORNL-6086/V1 (Oak Ridge, Tenn.: Oak Ridge National 
Laboratories) 

, 129 

Bertaux, 1. L., Kyrola, E., Quemerais, E., et al. 1995, Sol. Phys., 162, 403 
Bzowski, M. 2001, in The Outer Heliosphere: The Next Frontiers, ed. K. Scherer, 

H. Fichtner, H. J. Fahr, & E. Marsch. COSPAR Colloquia Series Vol. 1 1 
(Elsevier, Pergamon), 69-72 

Bzowski, M. 2003, A&A, 408, 1155 

Bzowski, M., Makinen, T., Kyrola, E., Summanen, T., & Quemerais, E. 2003, 
A&A, 408, 1165 

Bzowski, M., Mobius, E., Tamopolski, S., Izmodenov, V, & Gloeckler, G. 2008, 

ArXiv e-prints, 0710.1480, submitted to A&A 
Bzowski, M. & Tamopolski, S. 2006, in American Institute of Physics 

Conference Series, Vol. 858, Physics of the Inner Heliosheath, ed. 

I. Heerikhuisen, V. Florinski, G. P. Zank, & N. V. Pogorelov, 251-256 
Cook, J. W., Brueckner, G. E., & van Hoosier, M. E. 1980, J. Geophys. Res., 85, 

2257 

Cook, J. W., Meier, R. R., Brueckner, G. E., & van Hoosier, M. E. 1981, A&A, 
97, 394 

Fahr, H., Fichtner, H., & Scherer, K. 2007, Rev. Geophys., 45, RG4003, 

doi: 10. 1029/2006RG000214 
Fahr, H.-l. & Scherer, K. 2004, ASTRA, 1, 3 
Gruntman, M. 1997, Review of Scientific Instruments, 68, 3617 
Gruntman, M., Roelof, E. C, Mitchell, D. G., et al. 2001, J. Geophys. Res., 106, 

15767 

Gruntman, M. A. 1990, Planet. Space Sci., 38, 1225 
Gruntman, M. A. 1992, Planet. Space Sci., 40, 439 

Heerikhuisen, 1., Pogorelov, N. V., Zank, G. P, & Florinski, V. 2007, Apl, 655, 
L53 

Hsieh, K. C, Shih, K. L., lokipii, J. R., & Grzedzielski, S. 1992, ApJ, 393, 756 
Issautier, K., Meyer-Vemet, N., Moncuquet, M., & Hoang, S. 1998, 

1. Geophys. Res., 103, 1969 
King, I. H. & Papitashvih, N. E. 2005, 1. Geophys. Res., 110, 2104 
Kyrola, E., Summanen, T., Makinen, T, et al. 1998, J. Geophys. Res., 103, 14523 
Lemaire, P., Emerich, C, Vial, I.-C, et al. 2005, Advances in Space Research, 

35, 384 

Lemaire, P L., Emerich, C, Vial, 1. C, et al. 2002, in ESA SP-508: From Solar 

Min to Max: Half a Solai- Cycle with SOHO, 219-222 
Lindsay, B. G. & Stebbings, R. F 2005, 1. Geophys. Res., 1 10, 12213 
Lotz, W. 1967, Z. Phys., 206, 205 

Maksimovic, M., Gary, S. P, & Skoug, R. M. 2000, 1. Geophys. Res., 105, 18337 
McComas, D., Allegrini, F, Bartolone, L., et al. 2004, in American Institute of 

Physics Conference Series, Vol. 719, Physics of the Outer Heliosphere, 162- 

181 

McComas, D., Allegrini, F, Bartolone, L., et al. 2005, in ESA Special 

Publication, Vol. 592, Solar Wind 11/SOHO 16, Connecting Sun and 

Heliosphere, 689-692 
McComas, D. 1., Allegrini, F, Bartolone, L., et al. 2006, in American Institute 

of Physics Conference Series, Vol. 858, Physics of the Inner Heliosheath, ed. 

1. Heerikhuisen, V. Florinski, G. P Zank, & N. V. Pogorelov, 241-250 
McComas, D. 1., Barraclough, B. L., Funsten, H. O., et al. 2000, 

1. Geophys. Res., 105, 10419 
McComas, D. 1., Elliot, H. A., Schwadron, N. A., et al. 2003, 

Geophys. Res. Lett., 30, 10.1029/2003GL017136 
McComas, D. 1., Funstein, H. O., Gosling, J. T., & Pryor, W. R. 1999, 

Geophys. Res. Lett., 26, 2701 
Ogawa, H. S., Wu, C. Y. R., Gangopadhyay, P, & ludge, D. L. 1995, 

1. Geophys. Res., 100, 3455 



16 



M.Bzowski: Modifications of H ENA between termination shock and Earth orbit 



Phaneuf, R. A., Janev, R. K., & Pindzola, M. S. 1987, Atomic data for fusion. 
Vol. ORNL-6086A'5 (Oak Ridge, Tenn.: Oak Ridge National Laboratories) 

Phillips, J. L., Bame, S. J., Barnes, A., et al. 1995, Geophys. Res. Lett., 22, 3301 

Pilipp, W. G., Muehlhaeuser, K.-H., Miggenrieder, H., Montgomery, M. D., & 
Rosenbauer, H. 1987, J. Geophys. Res., 92, 1075 

Pryor, W. R., Ajello, J. M., Barth, C. A., et al. 1992, ApJ, 394, 363 

Quemerais, E. 2000, A&A, 480, 353 

Quemerais, E., LaUement, R., Koutroumpa, D., & Lamy, P. 2007, ApJ, 667, 1229 
Rucihski, D., Cummings, A. C, Gloeckler, G., et al. 1996, Space Sci. Rev., 78, 

73 

Rucihski, D. & Fahr, H. J. 1989, A&A, 224, 290 
Scherer, H. & Fahr, H. J. 1996, A&A, 309, 957 
Scherer, K. & Fahr, H.-J. 2003, A&A, 404, L47 

Scime, E. E., Bame, S. J., Feldman, W. C, Gary, S. P, & Phillips, J. L. 1994, 

J. Geophys. Res., 99, 23401 
Sternal, O., Fichtner, H., & Scherer, K. 2008, A&A, 477, 365 
Tamopolski, S. & Bzowski, M. 2007, ArXiv Astrophysics e-prints 0701133, sub- 
mitted to A&A 

Tobiska, W. K., Woods, T., Eparvier, R, et al. 2000, J. Atm. Terr. Phys., 62, 1233 
von Steiger, R., Schwadron, N. A., Fisk, L. A., et al. 2000, J. Geophys. Res., 

105, 27217 
Witte, M. 2004, A&A, 426, 835 

Woods, T. N., Tobiska, W. K., Rottman, G. J., & Worden, J. R. 2000, 
J. Geophys. Res., 105, 27195 



